Luteolin Causes 5′CpG Demethylation of the Promoters of TSGs and Modulates the Aberrant Histone Modifications, Restoring the Expression of TSGs in Human Cancer Cells

Cancer progression is linked to abnormal epigenetic alterations such as DNA methylation and histone modifications. Since epigenetic alterations, unlike genetic changes, are heritable and reversible, they have been considered as interesting targets for cancer prevention and therapy by dietary compounds such as luteolin. In this study, epigenetic modulatory behaviour of luteolin was analysed on HeLa cells. Various assays including colony forming and migration assays, followed by biochemical assays of epigenetic enzymes including DNA methyltransferase, histone methyl transferase, histone acetyl transferase, and histone deacetylases assays were performed. Furthermore, global DNA methylation and methylation-specific PCR for examining the methylation status of CpG promoters of various tumour suppressor genes (TSGs) and the expression of these TSGs at transcript and protein level were performed. It was observed that luteolin inhibited migration and colony formation in HeLa cells. It also modulated DNA methylation at promoters of TSGs and the enzymatic activity of DNMT, HDAC, HMT, and HAT and reduced the global DNA methylation. Decrease in methylation resulted in the reactivation of silenced tumour suppressor genes including FHIT, DAPK1, PTEN, CDH1, SOCS1, TIMPS, VHL, TP53, TP73, etc. Hence, luteolin-targeted epigenetic alterations provide a promising approach for cancer prevention and intervention.


Introduction
Cancer chemoprevention is the use of compounds that are natural, synthetic, or chemical to reverse, retard, or prevent the onset of cancer using mechanisms that kill malignant cells selectively and spare normal cells [1,2]. These cancer chemopreventive agents can act at all three stages of cancer development-namely, initiation, promotion, and progression. FDA approved chemopreventives such as tamoxifen and raloxifene have the major side effect of causing blood clotting; similarly, tamoxifen-if given for prolonged period-can lead to endometrial cancer. Hence, researchers are shifting their focus towards plant-derived chemopreventive agents, as they have few or no side effects, are more affordable, easily available, and well tolerated. Studies indicate that dietary bioactive compounds such as polyphenols, isothiocyanates, etc. are found to have a profound effect on cancer cells and act via targeting various hallmarks of cancer, leading to the reversal of carcinogenesis [1,3,4]. Various phenolic agents such as EGCG (epigallocatechin gallate), curcumin, etc. have been reported to exert antiproliferation, apoptosis-inducing, and antioxidant properties and to restore the function of TSGs by targeting various molecular targets and pathways [1,[4][5][6].
Cumulative data, both experimental and epidemiological, approves the utilization of dietary phytochemicals in chemoprevention, while suggesting that regular consumption of these phytochemicals is a promising new method for cancer prevention [7]. Evidence from many scientific reports have indicated that various bioactive compounds found in vegetables and fruits influence epigenetic modulations by regulating crucial processes that entail tumour suppressor genes reactivation, oncogene suppression, cell-survival protein activation, and apoptosis induction in various cancers [8]. Several nutrients in our daily dietary intake play an important role in DNA methylation either by modulating (DNA methyl transferase) DNMT expression or by altering the presence of methyl donors [8].
Dietary bioactive compounds such as curcumin, tea polyphenols, sulforaphane (SFN), genistein, resveratrol, and others have been found to be effective against cancer cells by altering regulatory mechanisms that directly affect the epigenome [9]. Bioactive agents act upon cancer progression and reverse the process by interfering with the action of (histone deacetylases) HDACs, DNMTs, and HATs (histone acetyl transferase) [8][9][10][11]. Published reports of polyphenol show that they are able to modify chromatin structure and reactivate silenced genes [12,13], potentially reverting malignancy-associated epigenetic alterations in cell lines [2,14]. Hence, modulation of epigenetic mechanisms poses a prominent target for cancer therapeutics [15]. This multidimensional molecular approach of polyphenols can be substantially highlighted for cancer prevention and therapy. Polyphenols have shown the inhibition of proliferation, invasion, and modulation of epigenetic pathways in different cell lines [16][17][18].
Luteolin, a flavone, has been studied as an anticancer agent and is known to inhibit tumour development in cancer cells by inducing apoptosis, inhibiting metastasis by modulation of different pathways such as AKT, MAPK, etc. [19][20][21][22][23][24]. Luteolin has also shown anti-invasive and epigenetic modulatory action such as modulation of methylation and inhibition of HDAC on different cell lines [17,25,26]. In addition to this, flavones have shown a remarkable safety profile (zero toxicity at up to 140 g/day) with no adverse side effects [21]. However, few studies are available on epigenetics modulation such as DNA methylation and histone modifications by DNMT, HMT, and HDAC after luteolin treatment [25,27,28]. Hence, this study was intended to provide experimental evidence of luteolin-mediated changes, targeting its antimigratory and chromatin modulatory nature on HeLa cells for cancer prevention and treatment.

Figure 1. (A)
. Luteolin treatment of HeLa cells at 10 and 20 µM for 48 h demonstrated profound decrease in percentage methylation in 5′ CpG promoter regions of TSGs, as compared with the untreated controls, in a dose-dependent manner ( * represents that data is statistically significant with p ≤ 0.05). (B). Luteolin decreased global DNA methylation in HeLa cells. The data are presented as mean of three independent experiments ± SD (* represents that data is statistically significant with p ≤ 0.05).(C). Luteolin decreased DNMT activity in HeLa cells in a concentration-dependent manner. The activity of the treated was compared with the untreated, and the values plotted are the mean of three experiments ± SD (* represents that data is statistically significant with p ≤ 0.05).

Luteolin Diminishes Global DNA Methylation of HeLa Cells
Treatment of HeLa cells with 5, 10, and 20 µM of luteolin for 48 h reduced the global methylation in HeLa cells in a dose-dependent manner. Luteolin decreased global DNA methylation to 57%, 50%, and 30% at 5, 10, and 20 µM treatment for 48 h, respectively, in comparison with control ( Figure 1B).

Luteolin Diminishes DNMT Activity in HeLa Cells
Incubation of the nuclear extract with luteolin decreased the DNMT activity in HeLa cells in a dose-dependent manner. Luteolin at 5, 10, and 20 µM reduced the DNMT activity by 35%, 57%, and 66%, respectively, in comparison with the untreated control ( Figure 1C).

Luteolin Modulates Migration/Inflammation Related and Tumour Suppressor Genes
The decreased methylation percentage of various tumour suppressor genes after luteolin treatment was further verified by checking the expression of these genes using quantitative RT PCR. It was observed that luteolin increased the expression of various TSGs and the relative quantity (RQ) values after 10 and 20 µM treatment of HeLa cells for 48 h, respectively, for VHL (RQ 7,9); TIMP3 (RQ 4.5, 5.6); TIMP4 (RQ 6, 10); PTEN (RQ 2, 4.0); RARB (RQ 1.78,3.5); SOCS1 (RQ 2.65, 3.6); FHIT (RQ 1.5, 2.8); and GSTP1 (RQ 1.5, 1.8) and at the same time decreased the expression of various migration related genes such as MMP2 (RQ 0.53, 0.31); MMP9 (RQ 0.34, 0.38); MMP14 (RQ 0.28, 0.07); SNAIL2 (RQ 0.56, CpG promoter regions of TSGs, as compared with the untreated controls, in a dose-dependent manner (* represents that data is statistically significant with p ≤ 0.05). (B). Luteolin decreased global DNA methylation in HeLa cells. The data are presented as mean of three independent experiments ± SD (* represents that data is statistically significant with p ≤ 0.05). (C). Luteolin decreased DNMT activity in HeLa cells in a concentration-dependent manner. The activity of the treated was compared with the untreated, and the values plotted are the mean of three experiments ± SD (* represents that data is statistically significant with p ≤ 0.05).

Luteolin Diminishes Global DNA Methylation of HeLa Cells
Treatment of HeLa cells with 5, 10, and 20 µM of luteolin for 48 h reduced the global methylation in HeLa cells in a dose-dependent manner. Luteolin decreased global DNA methylation to 57%, 50%, and 30% at 5, 10, and 20 µM treatment for 48 h, respectively, in comparison with control ( Figure 1B).

Luteolin Diminishes DNMT Activity in HeLa Cells
Incubation of the nuclear extract with luteolin decreased the DNMT activity in HeLa cells in a dose-dependent manner. Luteolin at 5, 10, and 20 µM reduced the DNMT activity by 35%, 57%, and 66%, respectively, in comparison with the untreated control ( Figure 1C).

Luteolin Modulates Migration/Inflammation Related and Tumour Suppressor Genes
The decreased methylation percentage of various tumour suppressor genes after luteolin treatment was further verified by checking the expression of these genes using quantitative RT PCR. It was observed that luteolin increased the expression of various TSGs and the relative quantity (RQ) values after 10 and 20 µM treatment of HeLa cells for 48 h, respectively, for VHL (RQ 7, 9); TIMP3 (RQ 4.5, 5.

Luteolin Decreases HDAC Activity in HeLa Cells
Histone deacetylation leads to heterochromatin and thus reduced gene expression, so the activity of histone deacetylases was analysed after direct incubation of nuclear extracts with different concentrations of luteolin. Incubation of the nuclear extract with 5, 10, and 20 µM luteolin led to decreased HDAC activity by 28%, 40%, and 55%, respectively, as compared with the controls, as shown in Figure 2B.

Luteolin Decreases HDAC Activity in HeLa Cells
Histone deacetylation leads to heterochromatin and thus reduced gene expression, so the activity of histone deacetylases was analysed after direct incubation of nuclear extracts with different concentrations of luteolin. Incubation of the nuclear extract with 5, 10, and 20 µM luteolin led to decreased HDAC activity by 28%, 40%, and 55%, respectively, as compared with the controls, as shown in Figure 2B.  Figure 2C). CIITA regulates the transcription of major histocompatibility complex (MHC) class I and II genes.

Luteolin Reduces the HAT Activity in a Dose-Dependent Manner
Luteolin at 5, 10, and 20 µM treatment decreased the HAT activity in HeLa cells, as compared with the untreated controls. On incubation of the nuclear extracts with different concentrations of luteolin, such as 5, 10, and 20 µM, there was a decrease in HAT activity by 17%, 51%, and 66%, respectively, in comparison with the untreated controls ( Figure 3A). Figure 2C). CIITA regulates the transcription of major histocompatibility complex (MHC) class I and II genes.

Luteolin Reduces the HAT Activity in a Dose-Dependent Manner
Luteolin at 5, 10, and 20 µM treatment decreased the HAT activity in HeLa cells, as compared with the untreated controls. On incubation of the nuclear extracts with different concentrations of luteolin, such as 5, 10, and 20 µM, there was a decrease in HAT activity by 17%, 51%, and 66%, respectively, in comparison with the untreated controls ( Figure  3A).

Luteolin Decreases the HMT H3K9 Enzyme Activity in HeLa Cells
Histone methyl transferase HMT H3K9 enzyme causes mono, di-, and tri-methylation at histone 3 lysine 9. These marks are repressors of transcription. H3K9 methylation and DNA methylation are cooperative epigenetic changes that complement each other and are responsible for silencing of the promoters of several genes. Luteolin treatment leads to inhibition of HMT by 24%, 31%, and 47% after incubation of the nuclear extract with 5, 10, and 20 µM ( Figure 3B), respectively. The inhibition percentage was obtained by comparing treated sample activity with the control sample activity.

Luteolin Decreases the HMT H3K9 Enzyme Activity in HeLa Cells
Histone methyl transferase HMT H3K9 enzyme causes mono, di-, and tri-methylation at histone 3 lysine 9. These marks are repressors of transcription. H3K9 methylation and DNA methylation are cooperative epigenetic changes that complement each other and are responsible for silencing of the promoters of several genes. Luteolin treatment leads to inhibition of HMT by 24%, 31%, and 47% after incubation of the nuclear extract with 5, 10, and 20 µM ( Figure 3B), respectively. The inhibition percentage was obtained by comparing treated sample activity with the control sample activity.

Luteolin Modulates Protein Expression of Genes Related to Migration, Inflammation, and TSGs
Proteome profiler-based quantitation of proteins involved in proliferation, migration, and various cellular events showed modulation after luteolin treatment, consistent with the mRNA expression. Luteolin resulted in downregulation of expression of various proteins related to migration and inflammation in HeLa cells, and fold changes (FC) after treatment with 10 and 20 µM for 48 h are given, respectively, for MMP2 (FC 0. 33 decreased after chrysin treatment ( Figure 3D). Phosphorylation marks of H3ser28p, H4ser1p, H4R3m2a, and H4Rm2s were also decreased after treatment of HeLa cells with luteolin 20 µM for 48 h ( Figure 3C,D).

Luteolin Modulates Protein Expression of Genes Related to Migration, Inflammation, and TSGs
Proteome profiler-based quantitation of proteins involved in proliferation, migration, and various cellular events showed modulation after luteolin treatment, consistent with the mRNA expression. Luteolin resulted in downregulation of expression of various proteins related to migration and inflammation in HeLa cells, and fold changes (FC) after treatment with 10 and 20 µM for 48 h are given, respectively, for MMP2 (FC 0. 33

Luteolin Repressed Colony Formation and Migration of HeLa Cells
Luteolin inhibits colony formation in HeLa cells. In untreated control, a plating efficiency of 94% was observed, and at 5 µM and 10 µM luteolin treatment for 48 h, 150 and 50 colonies were seen, and at 20 µM, almost no colonies were formed ( Figure 5A). However, the number of cells in untreated colonies was almost 10 to 20 cells per colony, whereas in treated colonies, the size was only 5-6 cells per colony, indicating that luteolin significantly hindered colony formation and induced the cytostatic behaviour in HeLa cells.

Luteolin Repressed Colony Formation and Migration of HeLa Cells
Luteolin inhibits colony formation in HeLa cells. In untreated control, a plating efficiency of 94% was observed, and at 5 µM and 10 µM luteolin treatment for 48 h, 150 and 50 colonies were seen, and at 20 µM, almost no colonies were formed ( Figure 5A). However, the number of cells in untreated colonies was almost 10 to 20 cells per colony, whereas in treated colonies, the size was only 5-6 cells per colony, indicating that luteolin significantly hindered colony formation and induced the cytostatic behaviour in HeLa cells.

Luteolin Inhibits Migration Capacity of HeLa Cells
The migration inhibitory effect of luteolin was further confirmed by transwell assay. The migration of control cells was used to compare the migration in treated samples. The width of the wound (in cm) gradually increased in the wells treated with 10 µM for 24 h and 48 h of luteolin by 7% and 23%, respectively, and increased by 15% and 31% after 20 µM treatment for 24 h and 48 h, respectively, while the control wells exhibited the inverse, with width decreasing in a time-dependent manner, by an average of 45% by 24 h and complete closure at 48 h ( Figure 5B). The migration inhibitory effect of luteolin was further confirmed by transwell assay. The migration of control cells was used to compare the migration in treated samples. Luteolin decreased the migration capacity of HeLa cells to 18% and 2% at 10 and 20 µM treatment for 48 h in comparison with the control ( Figure 5C).

Luteolin Inhibits Migration Capacity of HeLa Cells
The migration inhibitory effect of luteolin was further confirmed by transwell assay. The migration of control cells was used to compare the migration in treated samples. The width of the wound (in cm) gradually increased in the wells treated with 10 µM for 24 h and 48 h of luteolin by 7% and 23%, respectively, and increased by 15% and 31% after 20 µM treatment for 24 h and 48 h, respectively, while the control wells exhibited the inverse, with width decreasing in a time-dependent manner, by an average of 45% by 24 h and complete closure at 48 h ( Figure 5B). The migration inhibitory effect of luteolin was further confirmed by transwell assay. The migration of control cells was used to compare the migration in treated samples. Luteolin decreased the migration capacity of HeLa cells to 18% and 2% at 10 and 20 µM treatment for 48 h in comparison with the control ( Figure 5C).

Discussion
Cancer is one of the main causes of mortality globally, and metastasis plays a major role in prognosis of the disease. The key enzymes involved in metastasis are the metalloproteases, and these can be potential targets for inhibition of metastasis [29]. Cancer, being a complex disease, is considered to be a manifestation of both genetic and epigenetic modulations, resulting in altered expression of genes related to proliferation, apoptosis, cell cycle, and migration [30,31] Epigenetic modulations in cancer cells include global DNA hypomethylation and hypermethylation of specific promoter regions of tumour suppressor genes (TSGs) and genes related to apoptosis and cell migration [30,32]. In spite of the fact that synthetic inhibitors showed promising results in alteration of these epigenetic modifications-owing to poor oral availability, toxicity, and lack of specificity-the ongoing search for compounds that are more specific and less toxic continues [32,33]. Interestingly, polyphenols including flavonoids have shown a number of biological activities

Discussion
Cancer is one of the main causes of mortality globally, and metastasis plays a major role in prognosis of the disease. The key enzymes involved in metastasis are the metalloproteases, and these can be potential targets for inhibition of metastasis [29]. Cancer, being a complex disease, is considered to be a manifestation of both genetic and epigenetic modulations, resulting in altered expression of genes related to proliferation, apoptosis, cell cycle, and migration [30,31] Epigenetic modulations in cancer cells include global DNA hypomethylation and hypermethylation of specific promoter regions of tumour suppressor genes (TSGs) and genes related to apoptosis and cell migration [30,32]. In spite of the fact that synthetic inhibitors showed promising results in alteration of these epigenetic modifications-owing to poor oral availability, toxicity, and lack of specificity-the ongoing search for compounds that are more specific and less toxic continues [32,33]. Interestingly, polyphenols including flavonoids have shown a number of biological activities such as antiviral, antioxidant, and pro-apoptotic effects both in vitro and in vivo [34,35]. Additionally, several flavonoids such as apigenin and luteolin have shown inhibitory effects on different epigenetic enzymes [32,36] and genes related to migration [20,37,38]. However not much is known about the effect of luteolin on cervical cell migration via modulation of the epigenetic pathway. Hence, this study was planned to analyse the modulation of epigenome after luteolin treatment and to study the effect of luteolin on migration, inflammation, and growth of HeLa cells.
Luteolin has displayed anticancer effects such as induction of apoptosis and induction of cell cycle in different cell lines [39,40]. It has shown induction of apoptosis via inhibition of E6/E7 viral proteins in cervical cancer cells [40]. Similar reports of induction of apoptosis and inhibition of various pathways in HeLa cells after luteolin treatment were visualised in our earlier study [24], wherein it was seen that luteolin caused differential cytotoxicity towards HeLa cells in comparison with the lymphocytes in the concentration range (1-40 µM). There were no significant effects of the used luteolin concentration on lymphocytes [24]. The cell viability percentages at 10 and 20 µM were 70% and 50%, respectively. These two sublethal doses (10 and 20 µM) were taken for further studies. Luteolin inducted cell cycle arrest and induced apoptosis by decreasing the expression of AKTs, cyclins, CDK2, Bcl-2, and MCL1, as well as increasing the expression of Bax, caspase-3, 8, and 9, and p21 [24]. However, different levels of cytotoxicity have been observed in different cell lines after luteolin treatment [25,26,41,42]. Luteolin has shown selective cytotoxicity towards hepatocytes of hepatocellular carcinoma (HCC) rats, as compared with normal hepatocytes [43]. Additionally, some reports have attributed anticancer effect of luteolin to being sensitive against drug resistant cancer cells [44].
Interestingly some reports have shown inhibition of metastasis in colorectal cancer (HT-29 and SW480 and downregulation of MMP 2, MMP 3, and MMP 9 via inhibition of the AKT/PI3K pathway in melanoma cells (A375 cells) after luteolin treatment [37,38]). Luteolin has depicted increased expression of MiR-26a, which is a regulator of EZH2, and at the same time, it has inhibited EZH2 and H3K27me3, leading to self-programmed death and cycle arrest and apoptosis in prostate cancer cells [45]. To understand the epigenetic modulation by luteolin, DNA methylation and histone modulation was assessed after luteolin treatment. DNA methylation of TSGs is associated with silencing of the genes including GSTP1, BRCA1, RARB, MGMT, RASSF1A, and CDKN2A/p16 [46][47][48][49][50]. In the present study, it was found that luteolin decreased the DNA methylation percentage in TSGs ( Figure 1A). These TSGs are found to be hypermethylated in cervical cancer, hence leading to loss of their functions [51]. The decrease in methylation at promoters of TSGs was further affirmed by reduction in DNMT activity by 35%, 57%, and 66% after 5, 10, and 20 µM luteolin treatment ( Figure 1C). Luteolin also showed a decrease in global DNA methylation to 57%, 50%, and 30% after 5, 10, and 20 µM treatment for 48 h ( Figure 1B). The decrease in methylation of TSGs led to the reactivation of the genes including TP53, GSTP1, FHIT, SOCS1, RARB, PTEN, TP73, TIMP3, TIMP4, and VHL (Figure 2A) at the transcript level and E-cadherin and Tp53 at the protein level. Similar results of decrease in DNMT activity, global DNA methylation, and reactivation of various TSGs has been demonstrated by other phytochemicals, including luteolin [26,52].
Gene expression is also modulated by histone modifications that lead to hetero chromatin or euchromatin formation. Histone modifications are mediated by epigenetic enzymes such as histone acetylases, deacetylases, methylases, and demethylases which are found to be aberrantly expressed in cancer cells [53]. Luteolin showed direct inhibition of the biochemical activity of HDAC, with a sharp decline of 28%, 40%, and 55% at concentrations of 5, 10, and 20 µM ( Figure 2B). This was well corelated with the downregulation of HDAC1, HDAC5, HDAC7, HDAC2, HDAC10, and HDAC11 at transcript level in a dose-dependent manner ( Figure 2C). The results by other scientists on human epithelioid cancer cells have demonstrated luteolin to be a potent HDAC inhibitor [25,30]. Further, it was seen that luteolin treatment decreased the expression of HAT 1 at the transcript level ( Figure 2C) and that HAT enzymatic activity also decreased progressively after luteolin treatment ( Figure 3A). Biochemical activity of HMT (H3K9) ( Figure 3B Figure 2C). These enzymes have been implicated in specific stages of cell migration, cell propagation, or metastasis, as described as follows. ESCO2 has the ability to knock down MMP2, thus inhibiting migration of cancer cells [54]. SETD2 is a tumour suppressor leading to trimethylation of H3K36me3, which is a repressive mark [55]. Major histocompatibility complex class II is linked to better prognosis of cancer patients, and its expression on cancer cells depends on the trans activator CIITA, which is inactivated in colorectal and gastric cancers [56], hence attributing tumour suppressor function to CIITA [53]. WHSC1 activates TWIST and thus EMT; hence, its downregulation inhibits migration [57]. Inhibition of Aurora kinases has been linked to suppression of cell propagation and metastasis [58,59] and to inhibition of the development and progress of many cancers [60,61]. ESCO1 expression is also linked to cancer progression and metastasis [62]. Hence, downregulation of ESCO1, WHSC1, AURKA, AURKB, and AURKC by luteolin points towards its antiproliferative and antimigratory role.
The modulation of the proliferative and migratory pathway was further verified by colony assay, scratch-wound, and invasion assays after 5, 10, and 20 µM treatment of luteolin. The number of colonies formed decreased from 480 in control to 150 and 50 at 5 and 10 µM for 48 h; however, no colonies could be seen at 20 µM treatment ( Figure 5A), indicating complete cytostatic induction. In scratch wound assay, the increase in wound width was visualized by 7% and 23% at 10 µM treatment for 24 and 48 h and 15% and 31% increase in wound width at 20 µM luteolin treatment for 24 and 48 h treatment, respectively ( Figure 5B). This was further supported by transwell assay, wherein luteolin decreased the migration capacity of HeLa cells to 18% and 2% at 10 and 20 µM treatment for 48 h ( Figure 5C). These findings further supported the downregulation of pro-migratory (MMPs) genes and upregulation of antimigratory (TIMPs and PTEN) and antiproliferative genes. Similar results of inhibition of migration and invasion by luteolin in colorectal cells and squamous cell carcinoma cells were obtained by some authors [37,38,73].
Our study systematically demonstrates the mode of action of luteolin, which involves downregulating the expression and activity of epigenetic modulators such as DNMTs, HDACs, HMTs, and HAT, leading to consequential demethylation of promoters of TSGs and reactivation of TSGs. In addition, luteolin downregulated the expression of genes related to proliferation and migration, such as MMPS, SNAIL, MTA1, and MTA2, and upregulated inhibitors of migration such as TIMPs and TSGs such as PTEN.

Maintenance of Cervical Cancer Calls (HeLa) and Drug Dilution
Human cervical cancer cell line (HeLa) was maintained in DMEM supplemented with 10% FBS and incubated at 37 • C. Luteolin was purchased from Sigma-Aldrich (Merck, KgaA, Darmstadt, Germany) and a stock of 69.87 mM was prepared with DMSO. Further concentration of 1mM was prepared and concentrations (1-40 µM) were made.

Methylation-Specific PCR (MSRE-PCR)
An EpiTect Methyl II PCR kit (catalogue number 335452, Qiagen, Germantown, MD, USA) was procured. This kit helps in the quantitation of methylated and unmethylated DNA. Restriction digestion was conducted on 1 µg of genomic DNA (treated and DMSOtreated), and the cleaved DNA in each reaction was measured using quantitative PCR. The amount of unmethylated and methylated DNA was assessed using the ∆∆ CT technique, which compared the quantity of each reaction to that of a mock (no enzymes added) response. The gene panel (which included pre-designed primers) included tumour suppressor genes. The mean of three experiments was used to assess statistical significance (p < 0.05).

Global DNA Methylation Assay
The DNA of the treated (5, 10, and 20 µM of luteolin for 48 h) and vehicle control cells was extracted using the GenElute Mammalian Genomic DNA Miniprep Kit (Catalogue number G1N70; Sigma-Aldrich, Merck, KgaA). The quality and quantity of DNA was checked by gel electrophoresis (1% agarose gel, Catalog No. A9539, Sigma-Aldrich, Merck, KgaA) and spectrophotometry using Nanodrop 2000 (Thermo Scientific™, Waltham, MA, USA), respectively. The MethylFlash™ Methylated DNA Quantification Kit (Catalogue number: P-1034; Epigentek, Farmingdale, NY, USA) was used for the detection of methylated DNA by 5-mC antibody, which can be assessed calorimetrically.
The treated and untreated DNA were investigated for global DNA methylation as per the protocol. Absorbance was read on an ELISA reader at 450 nm, 10 min after placement in developing solution. Optical density values are relative to the quantity of methylated DNA, and the levels of methylation were calculated in comparison with the control. Data were taken as the average of 3 experiments ± SD. Significant differences were established at p ≤ 0.05.

Expression Analysis of TSGs
GenElute Mammalian Genomic Total RNA Kit, Sigma, Merck, KgaA, was procured. The manufacturer's instructions were followed while extracting RNA from luteolin-treated HeLa (10 and 20 M for 48 h) and untreated cells. The Applied Biosystems High-Capacity cDNA Reverse Transcription Kit was then used to synthesize cDNA on the RNA. The expression of various TSGs and genes encoding various pathways was analysed using a TaqMan-based custom array (4391524). Thermo Fisher's DataAssistTM software was used to perform the PCR array on QuantStudio3 and analyse it using the ∆∆ CT technique. ACTB (Actin Beta) expression was used for normalization.

Protein Expression by Proteome Profiler Array
Proteome profiler array (catalogue no. ARY026) was procured from R&D, USA. The treated (10 and 20 µM for 48 h) and untreated cells were collected and resuspended in lysis buffer 17 (1 mL per 10 7 cells); the lysate was prepared, and protein concentration was checked using Pierce BCA Assay (Catalogue no: 23225; Thermo-Fisher Scientific, Waltham, MA, USA). The signals produced on the membranes were quantified by chemiluminescent gel doc system (Bio-Rad Laboratories, Hercules, CA, USA) and analysed using Image Lab software (version 6.1). Significance was calculated by one-way ANOVA using SPSS software and was documented at p ≤ 0.05.

Analysis of Epigenetic Enzymes Involved in Chromatin Modification
Human Chromatin Modification Enzymes Array (Catalogue number PAHS-085Z; Qiagen, MD, USA) was used to check the expression of the genes such as DNMTs, HMTs, HATS, and HDMs. RNA at a concentration of 1 µg was used to synthesize cDNA, and it was diluted to 1350 µL with nuclease free water, and an equal amount of RT2 SYBR ® Green qPCR Master mix (Catalogue number: 330504; Qiagen, MD, USA) was added to this. From this mixture, 25 µL was poured into each well of an array plate having pre-defined primers, and the plate was run on ABI Quant Studio 3. The normalization was performed using GAPDH as the housekeeping gene, and fold change was calculated by comparing treated samples with untreated controls. The statistical significance was calculated, and the significance was established at p ≤ 0.5.

Nuclear Extract Preparation
Nuclear extracts from HeLa cells were prepared using EpiQuik TM nuclear extraction kit (Catalogue number OP-0002; Epigentek, NY, USA) as per the manufacturer's protocol. Approximately 2 × 10 6 untreated cells were trypsinised, collected, and centrifuged at 150× g for 10 min. The nuclear extract was prepared and stored at −80 • C until used. The protein concentration of the nuclear extract was checked using the BCA protocol.

DNMT Activity Assay
Nuclear extract from untreated HeLa cells was prepared using EpiQuik TM nuclear extraction kit (Catalogue number OP-0002; Epigentek, NY, USA). The EpiQuik DNMT activity assay kit (Catalogue number P-3001; Epigentek, NY, USA) was used to check the effect of the luteolin on DNMT activity. In the sample wells AdoMet (methyl group donor), nuclear extracts, and different concentrations of luteolin (5, 10, and 20 µM in assay buffer) were used to check for DNMT activity by following the protocol, whereas assay buffer was used in controls instead of the inhibitor. After incubating the plate for 2 h and developing the signal, the plate was read at 450 nm. The percentage of inhibition in comparison with control was calculated and plotted as a graph. The experiment was repeated thrice and the mean ± SD was used to plot a graph. One-way ANOVA was used to check for statistical significance, and p was established at ≤0.05.

HDAC Activity Assay
The EpiQuik HDAC activity assay kit (Catalogue number P-4002; Epigentek, NY, USA) was used to determine the effect of luteolin treatment on HDAC activity. Luteolin (5 µM, 10, and 20 µM), nuclear extract, and HDAC assay buffer were added to the assay plates and incubated for one hour at 37 • C to allow the action of the enzyme, followed by incubation with capture and detection antibody. In the controls, luteolin was substituted by assay buffer. The plate was read on an ELISA reader at 450 nm. Inhibition of HDAC activity by luteolin was calculated in comparison with control wells. Experiments were repeated thrice, and the mean ± SD was used to plot as a graph (p ≤ 0.05).

HMT H3K9 Activity Assay
The EpiQuik HMT H3K9 activity assay kit (Epigentek, NY, USA) was used to check the effect of luteolin on HMT H3K9 activity by following the manufacturer's protocol. Briefly, to histone 3 lysine substrate-coated assay plate biotinylated substrate, assay buffer and AdoMet were added and incubated for 1.5 h at 37 • C. In the case of test samples, luteolin at 5, 10, and 20 µM was added to the wells. After developing the colour, the absorbance was read on an ELISA reader at 450 nm. The inhibition of HMT by luteolin in comparison with the controls was calculated and plotted as a graph. The results were taken as the mean of 3 experiments ± SD (p ≤ 0.05).

HAT Activity
HAT activity/inhibition kit (catalogue number P-4003, Epigentek, NY, USA) was procured. Briefly, the nuclear extracts were incubated with the substrate for HAT for 60 min, followed by washing. Luteolin concentrations were added in the test samples.
The experiment was carried out as per the manufacturer's protocol. The plate was read at 450 nm on an ELISA plate reader, and the HAT activity inhibition was calculated from the average of three experiments. The results were calculated as the mean of three experiments ± SD, (p ≤ 0.05).

Scratch-Wound Assay
To understand the antimigratory effect of luteolin, the scratch-wound assay was carried out. The protocol was adapted from our lab [74]. Briefly, almost 5 × 10 5 cells were plated in a six-well plate and incubated at 37 • C overnight. Next day a "wound" or a cellfree line was created on a confluent monolayer of cells by scratching the monolayer with a pipette tip, and the cells were incubated in the presence of different luteolin concentrations (0, 10, and 20 µM), and the "healing" of the wound, which occurs through cell migration and growth towards the cell free zone, was monitored after 24 h and 48 h. An inverted microscope was used to capture images of the wounds in each well prior to treatment and after 24 h and 48 h of the treatment. The wound width was measured using the paint software. The change in wound width in comparison with the control was calculated. The experiment was triplicated, and results expressed ± SD (p < 0.05).

Colony Forming Assay
Approximately 2.5 × 10 5 cells were plated in a six well plate and followed by treatment with different concentrations of luteolin (5, 10, and 20 µM) for 48 h. Treated cells were collected after 48 h treatment, counted, and plated at approximately 500 cells/well. The cells were allowed to grow for two weeks. Then, cells were fixed in 100% methanol and stained with 0.1% crystal violet. The stained colonies were counted and images were taken using an Olympus inverted microscope (Labomed, Los Angeles, CA, USA), and ImageJ software was used to count the colonies [75].

Trans Well Chamber Assay
Briefly, 5.0 × 10 3 cells/well (both treated and untreated) were seeded on the upper side of the insert (8 µm pore size, BD Biosciences, Franklin Lakes, NJ, USA), and on the lower side of the insert in the 24-well plate, complete media was kept. Cells on the inner side of the insert were stained with 0.1% crystal violet, and cells on the outer side of the insert were scraped. Invasive potential of treated and untreated cells was compared, and the migrated cells were counted by using an inverted microscope (Olympus Corporation, Tokyo, Japan, 200× magnification). Graphs were plotted after repetition of the experiment three times, and the results are expressed as mean ± SD (p < 0.05).

Statistical Analysis
The results of the experiments are expressed as average ± standard deviation of three experiments. SPSS software version 21.0 was used to check the statistical significance with one-way ANOVA followed by Tukey's HSD post hoc analysis.

Conclusions
Luteolin acts as a potent antiproliferative, antimigratory agent. It decreased the 5 CpG promoter methylation of various TSGs, leading to re-expression of these TSGs and inactivated genes related to migration and inflammation. Considering its antimigratory and epigenetics modulatory behaviours warrants its use as a potential chemopreventive agent for therapeutic use, but only after clinical trials.

Data Availability Statement:
The datasets used and/or analysed during the present study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.